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Abstract. We highlight and discuss the importance of accounting for nebular emission in the 
SEDs of high redshift galaxies, as lines and continuum emission can contribute significantly 
or subtly to broad-band photometry. Physical parameters such as the galaxy age, mass, star- 
formation rate, dust attenuation and others inferred from SED fits can be affected to different 
extent by the treatment of nebular emission. 

We analyse a large sample of Lyman break galaxies from z ~ 3-6, and show some main 
results illustrating e.g. the importance of nebular emission for determinations of the mass-SFR 
relation, attenuation and age. We suggest that a fairly large scatter in such relations could be 
intrinsic. We find that the majority of objects (~ 60-70%) is better fit with SEDs accounting 
for nebular emission; the remaining galaxies are found to show relatively weak or no emission 
lines. Our modeling, and supporting empirical evidence, suggests the existence of two categories 
of galaxies, "starbursts" and "post-starbursts" (lower SFR and older galaxies) among the LBG 
population, and relatively short star- format ion timescales. 
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1. Introduction 

Until recently the importance of nebular emission on the SED of high redshift galaxies 
and hence on on our understanding of the physical parameters of these objects was 
overlooked. For example, when Spitzer started to detect the rest-frame optical part of 
the spectrum of high-redshift galaxies redshifted into broad-band filters of its IRAC 
camera at 3.6 and 4.5 /im, it was soon noted that some of these galaxies showed fairly 
large Balmer breaks (e.g. Egami et al. 2005, Eyles et al. 2005). Several authors interpreted 
the large break as being due to the presence of "old", evolved stellar populations with 
ages of several hundred Myr (Eyles et al. 2005, Labbe et al. 2010). If true this would 
imply the presence of very early galaxy populations, and very large formation redshift s 
(z/~10-30). 

However, since emission lines are ubiquituous in star- forming galaxies and since atomic 
physics causes most of the strong observed lines to be in the (rest-frame) optical part 
of the spectrum, longward of the Balmer break, nebular emission affects broad-band 
photometry differentially, and can cause an apparent Balmer break. The importance of 
this effect for high-z galaxies was pointed out by Schaerer & de Barros (2009), who 
showed that significantly younger ages could be obtained when accounting for nebular 
lines. 

It has now become clear (Schaerer & de Barros 2010, Ono et al. 2010, Lidman et al. 
2011) that nebular emission (both lines and continuum emission) must be taken into 
account for the interpretation of photometric measurements of the SEDs of star-forming 
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galaxies such as Lyman-alpha emitters (LAE) and Lyman break galaxies (LBGs) - the 
dominant galaxy populations at high-z. 

Furthermore, as testified by the presence of Lya emission, a large (and possibly grow- 
ing) fraction of the currently know population of star-forming galaxies at high redshift 
shows emission lines (at least Lya!; cf.Ouchi et al. 2008, Stark et al. 2010 and others). 

In parallel a lot of diverse evidence for galaxies with strong emission lines and/or 
strong contributions of nebular emission to broad-band filters has been found at different 
redshifts, e.g. by Shim et al. (2011), McLinden et al. (2011), Atek et al. (2011), Trump 
et al. (2011), van der Wei et al. (2011). 

Here we briefly summarise our method to account for nebular emission in a SED plus 
photometric-redshift fitting code and present some results obtained from a systematic 
study of a large sample of z ~ 3-6 LBGs. 



2. SED modeling 

Several evolutionary synthesis models, such as PEGASE, GALEV and others include 
nebular emission (cf. Fioc et al. 1999, Chariot & Longhetti 2001, Anders & Fritze 2003, 
Zackrisson et al. 2008). However, they were not used to analyse observations of high 
redshift galaxies. 

Our SED fitting tool, described in Schaerer & de Barros (2009, 2010), is based on 
a version of the Hyperz photometric redshift code of Bolzonella et al. (2000), modified 
to taking into account nebular emission. In de Barros et al. (2011, hereafter dBll) we 
consider a large set of spectral templates based on the GALAXEV synthesis models of 
Bruzual & Chariot (2003), covering different metallicities and a wide range of star forma- 
tion histories. A Salpeter IMF is adopted. Nebular emission from continuum processes 
and numerous emission lines is added to the spectra predicted from the GALAXEV mod- 
els as described in Schaerer & de Barros (2009), proportionally to the Lyman continuum 
photon production. The intergalactic medium (IGM) is treated with the prescription of 
Madau (1995). 

The free parameters of our SED fits are: redshift z, metallicity Z (of stars and gas), 
star formation (SF) history described by the timescale r (i.e. SFR ex exp(— t/r)), the 
age t* defined since the onset of star-formation, and attenuation Ay described by the 
Calzetti law. In some models we also vary the strength of the Lya line from zero to its 
maximum predicted by case B (see Schaerer et al. 2011). In dBll we consider 3 sets 
of SF histories, exponentially declining SF with variable timescales, SFR=const, and 
rising star- formation following Finlator et al. (2011). We use Monte-Carlo simulations to 
determine probability distribution functions for each parameter /quantity. 



3. Impact of nebular emission on physical parameters of star-forming 
galaxies from z ~ 3 to 8 

In Schaerer & de Barros (2009) we have demonstrated the effect nebular emission 
may have on age determinations of high-z galaxies. In Schaerer & de Barros (2010) we 
have applied our SED fitting tool to z ~ 6-8 galaxies recently observed with the new 
WFC3 camera onboard HST and other, more luminous galaxies. We have shown an 
overall trend of the star-formation rate (SFR) increasing with stellar mass, albeit with a 
large scatter. Also, we have shown that the specific SFR may be higher than previously 
thought, depending on model assumptions. We have also found indications for dust in 
some galaxies at this redshift, and a possible trend of increasing dust attenuation with 
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galaxy mass. However, since only few of these galaxies are detected beyond ^ 1.6-2 /im, 
the uncertainties of the physical parameters are fairly large. 

To improve on this aspect, and to examine systematically the effects of nebular emission 
on the physical parameters of distant galaxies, we have examined a large sample of z ~ 
3-6 galaxies. Some results from this study (dBll) are highlighted here. 

3.1. z ~ 3-6 galaxies: photometric data and selection 
We have used the GOODS-MUSIC catalogue of Santini et al. (2009) providing photome- 
try in the U, B435, V606? 1776? Z 850LP, J, H, K, bands mostly from the VLT and HST, and 
the 3.6, 4.5, 5.8, and 8.0 /im bands from the IRAC camera onboard Spitzer. Using stan- 
dard criteria we have selected U, B, V, and i-drop galaxies. To reduce the contamination 
rate (typically ~ 10-20 %) we have only retained the objects whose median photometric 
redshifts agree with the targetted redshift range. This leaves us with a sample of 389, 
705, 199, and 60 galaxies with median photometric redshifts of %>hot = 3.3, 3.9, 4.9, and 
6.0. See dBll for more details. 

3.2. Results - Two categories of high-z galaxies 

One of the first results obtained when applying SED models with/without nebular emis- 
sion to our large galaxy sample is that the majority of the objects (~ 60-70 %) are 
better fit with nebular emission, irrespectively of the SF history. For the rest, the dif- 
ference between fits with or without nebular emission is generally minor. The difference 
between these two categories corresponds overall to objects whose best-fit SED clearly 
shows nebular lines (for the majority), and weak or no lines for the rest. Empirically this 
can be verified for z ~ 3.8-5 galaxies, where Ha falls in the 3.6 /im filter, whereas the 4.5 
/im filter remains essentially free of strong emission lines. Indeed, for this redshift range 
we find that the objects best fit with nebular lines show a systematically bluer (3.6-4.5) 
color than the other category (see de Barros et al. 2011b). Our first category strongly 
resembles the "Ha emitters" identified by Shim et al. (2011) through their 3.6 jam excess. 

From our modeling we suggest that the two categories correspond to "starburst" and 
"post-starburst" phases. This is supported by their difference in (instantaneous) SFR, 
where the objects better fit with emission lines ( "starbursts" ) show on average a higher 
SFR than the "post-starburst" category (see Fig. [1]), and also by an age difference. 
Interestingly we find a quite similar distribution of galaxies between these categories for 
all redshifts between z ~ 3 and 6. 

3.3. Results - Mass-SFR relations and the specific SFR at high-z 

The mass-SFR relation of galaxies and its evolution with redshift has drawn of lot of 
attention during the last few years. Some of the issues debated include the scatter in this 
relation and its behaviour at z ^> 3 (e.g. Daddi et al. 2007, Bouche et al. 2010). 

In Fig. [1] we show the derived SFR from our z ~ 4 sample as a function of UV 
magnitude and stellar mass respectively, and how these relations depend on assumptions 
on the star-formation histories and on the treatment of nebular emission. The reference 
model (left panels) assumes SFR=const and ages ^ 50 Myr, as often assumed in other 
studies. The derived SFR values are found on and above the SFR(UV) calibration of 
Kennicutt (1998) indicated by the dashed line, deviations from this line being due to 
dust attenuation. When nebular emission is treated and the SF history is kept free 
(among the declining SF histories) both effects contribute to increasing the scatter around 
the Kennicutt relation (right panels). Models with nebular emission naturally reveal 
a difference in SFR between the two categories discussed above, consistent with the 
observed color in the IRAC bands (3.6-4.5 jam). 
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Figure 1. SFR versus absolute UV magnitude (top) and SFR versus stellar mass (bottom) for 
the B-drop sample (z ~ 4) obtained from models with constant SFR without nebular emission 
(left panels) and for declining SFH and nebular emission (right panels). Points show the median 
values for each galaxy, color grades shows the 2D probability distribution function, the white 
solid line encloses the 68 % confidence limit. White points show the "starburst" objects, black 
"post-starburst" galaxies. The dashed lines show the SFR(UV) relation from Kennicutt (1998) 
and the mean SFR-M* relation of Daddi et al. (2007). 



The derived stellar masses do not strongly depend on the treatment of nebular emission. 
On average we obtain a reduction of ~ 30-40 % for z ~ 3-4 and somewhat more (~ 70-80 
%) at higher redshift. 

As shown in Fig. [1] the existence of a tight mass-SFR relation relies on the assumption 
of constant SFR over long enough timescales (^> 100 Myr). Again for variable SF histories 
and including nebular emission we find a large scatter (e.g. around the mean relation 
obtained by Daddi et al. (2007) at z ~ 2 shown in this Fig.). The median specific SFR 
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(=SFR/M* ) of z ~ 3-5 galaxies is typically a factor ~ 2-3 higher when nebular emission 
is accounted for both for SFR=const and for variable SFHs (cf. de Barros et al. 2011b). 

3.4. Results - Dust attenuation, galaxy ages and SF histories 

Our analysis also yields new information on dust attenuation and on trends with UV 
magnitude, galaxy mass, redshift etc. For illustration, we show the attenuation Ay versus 
^1500 for the U-, B-, and i-drop samples in Fig. [2j A clear evolution of decreasing Ay 
with redshift is found, both at a given UV magnitude and for the sample average. The 
exact value of the attenuation depends on model assumptions such as the SF history and 
of course also on the attenuation law. Taking nebular emission into account does not lead 
to a systematic shift of the inferred attenuation. 

At z rsj 3-4 the median age of galaxies shows a wide range from very young to close to 
the maximum age allowed by cosmology. Including the effects of nebular emission yields 
on average somewhat younger ages (as already found in our earlier papers), but does not 
preclude old ages. No strong evolution of the average galaxy age with redshift is found 
between z ~ 3 and 5 e.g. at fixed UV magnitude or given stellar mass; at higher z our 
sample appear significantly younger, though. 

Interestingly declining SF histories are clearly favoured for the majority of galaxies, 
as they fit better than models assuming SFR=const or rising SF. As already mentioned 
above, declining SF histories also provide a natural explanation for the two galaxy cat- 
egories identified ("starburst" and "post-starburst"). Other arguments, e.g. from galaxy 
statistics (Lee et al. 2009, 2011) or others (Stark et al. 2009), also tend to favour short 
times scales of SF. Timescales and SF histories of high-z galaxies are, however, subject 
to debate (e.g. Maraston et al. 2010, Finlator et al. 2011). For more detailed results and 
discussion of other implications see dBll . 




Figure 2. Attenuation Ay as a function of UV magnitude for the z ~ 3 (left), 4 (middle), and 
6 (right) samples, as determined from models with declining SFHs and accounting for nebular 
emission. Symbols as in Fig. 1. Note the fairly large scatter of Ay at given Muv, and the clear 
decrease of the attenuation at high redshift. 



4. Conclusions 

As our previous results on 2; ~ 6-8 galaxies (Schaerer& de Barros 2009, 2010), and 
the latest study of z ~ 3-6 objects (dBll) show, it is important to include the effects of 
nebular emission in SEDs models when deriving physical parameters from broad-band 
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photometric data. SED models including nebular lines and variable Lya can also be used 
to infer trends of Lya emission from photometric data, as shown in Schaerer et al. (2011. 

Accounting for nebular emission can lead to significant differences in the inferred age 
(which tend to become younger), masses (slightly lower), SFR , specific star- formation 
rate, attenuation and other parameters. We also find that the physical parameters depend 
quite significantly on the assumed star- formation histories. Our models favour relatively 
short (or recurrent?) star-formation timescales. 

From our analysis of a large sample of high-z galaxies we find that the majority of 
objects (~ 60-70%) is better fit with SEDs accounting for nebular emission (with a 
contribution proportionaly to the Lyman continuum flux). This category of apparently 
strongly SF objects ( "starburst" ) can also be distinguished empirically from the rest of 
the LBG population, which are probably closer to "post-starbursts" . We suggest that a 
fairly large scatter could be intrinsic in relations between SFR-mass and others. More 
detailed results and implications are presented and discussed in dBll. 
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